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ABSTRACT 

Relativistic outflows in the form of jets are common in many astrophysical objects. 
By their very nature, jets have angle dependent velocity profiles, T — F(r, 9, (ft), where 
T is the outflow Lorentz factor. In this work we consider photospheric emission from 
non-dissipative jets with various Lorentz factor profiles, of the approximate form F m 
ho/[($/$j) p + 1]j were 9- 3 is the characteristic jet opening angle. In collimated jets, 
the observed spectrum depends on the viewing angle, 9 V . We show that for narrow 
jets (^jTo < few), the obtained low energy photon index is a ss — 1 (dN/dE oc E a ), 
independent of viewing angle, and weakly dependent on the Lorentz factor gradient 
(p). A similar result is obtained for wider jets observed at 9 V s=s 9 y This result is 
surprisingly similar to the average low energy photon index seen in gamma- ray bursts. 
For wide jets (9jT > few) observed at 9 V -C 6j, a multicolor blackbody spectrum 
is obtained. We discuss the consequences of this theory on our understanding of the 
prompt emission in gamma-ray bursts. 

Key words: gamma-rays: bursts — plasmas — radiation mechanisms: thermal - 
radiative transfer — scattering — relativity 



1 INTRODUCTION 

Photospheric emission from highly relativistic outflows was 
early considered as an explanation for prompt gamma-ray 
bursts (GRBs, Goodman 1986; Paczyriski 1986). It is a nat- 
ural consequence of the fireball model, where the optical 
depth at the base of the outflow is much larger than unity 
(e.g. Piran 1999). Moreover, photospheric emission provides 
a natural explanation to the small dispersion of the sub-MeV 
peak and to the high prompt emission efficiency observed 
(Meszaros & Rees 2000). However, the observed spectrum 
usually appears significantly broader than a Planck spec- 
trum (Preece et al. 2000), being well fitted by a smoothly 
broken power law (the Band function, Band et al. 1993). 
Thus, the prompt emission has commonly been associated 
with synchrotron emission originating from kinetic energy 
dissipation outside the photosphere (Rees & Meszaros 1994). 
However, in recent years it has become clear that optically 
thin synchrotron emission is incompatible with the hard low 
energy slopes observed in a substantial fraction of GRBs 
(Preece et al. 1998; Kaneko et al. 2006; Beloborodov 2012). 
This has raised the need for alternative ideas. 

One appealing idea is broadening of the thermal spec- 
trum emitted from the photosphere. The emerging photon 
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spectrum from a static, optically thick, relativistic outflow 
can be widened in two ways: 

(i) Energy dissipation below the photosphere can heat 
electrons above the equillibrium temperature. These elec- 
trons emit synchrotron emission and Comptonize the ther- 
mal photons, thereby modifying the Planck spectrum (Rees 
& Meszaros 2005; Pe'er et al. 2005, 2006). The dissipation 
can be caused by shocks (Rees & Meszaros 2005; Lazzati 
et al. 2009; Ryde et al. 2011), dissipation of magnetic energy 
(Thompson 1994; Spruit et al. 2001; Giannios & Spruit 2005; 
Zhang & Yan 2011) or collisional processes (Beloborodov 
2010; Vurm et al. 2011). 

(ii) The photospheric radius is angle dependent 
(Abramowicz et al. 1991; Pe'er 2008). Moreover, it 
was shown by Pe'er (2008) that photons make their last 
scatterings at a distribution of radii and angles. The 
observer sees simultaneously photons emitted from a large 
range of radii and angles. Therefore, the observed spectrum 
is a superposition of comoving spectra. The Doppler boost 
is a function of angle, and the comoving temperature 
decreases with radius through adiabatic cooling. Depending 
on the outflow properties, this geometrical broadening can 
form observed spectra which appears significantly different 
from the Planck spectrum. 

Photospheric emission in the context of spherically sym- 
metric outflows has been studied by several authors. Good- 
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man (1986) considered a highly relativistic outflow in the 
context of cosmological GRBs. It was realized that the ob- 
served spectrum is broader than blackbody, however the 
analysis was one-dimensional. Abramowicz et al. (1991) re- 
alized that the two-dimensional shape of the photosphere 
in a relativistic, spherically symmetric wind is in fact con- 
cave and symmetric around the line-of-sight (LOS). This 
can be understood as a consequence of the dependence on 
viewing angle of the optical depth of relativistically moving 
media. Pe'er (2008) found a simple expression for the pho- 
tospheric radius, R p h oc (8 2 /3 + 1/r 2 ), where is the angle 
measured from the LOS, F = (1 — /3) -1 / 2 is the outflow 
bulk Lorentz factor and /3 is the outflow speed in units of 
the speed of light. Pe'er (2008) extended the photospheric 
emission model by recognizing the importance of considering 
photons from the entire emitting volume, introducing proba- 
bility density distributions for the last scattering photon po- 
sitions. Beloborodov (2011) took a different approach, solv- 
ing the radiative transfer equation in the relativistic limit. 
The approximate probability densities used by Pe'er (2008) 
was validated by Beloborodov (2011). 

All above mentioned works considered spherical explo- 
sions. As we show below, this is a good approximation for 
collimated outflows as long as the characteristic jet opening 
angle is much larger than 1/r and the outflow is observed 
at viewing angles much smaller than the jet opening angle. 
Within the collapsar model (MacFadyen & Woosley 1999) 
the jet is collimated by the pressure of the surrounding gas 
as it drills its way through the collapsing progenitor star. 
In such jets the position of the observer relative to the jet 
axis can affect the observed spectrum. Here we develop the 
theory of photospheric emission in collimated outflows, cal- 
culating the expression for the observed spectrum at any 
viewing angle. 

The mechanism responsible for jet collimation is not 
fully understood. Hydrodynamic simulations of GRB jets 
after the launching phase (e.g., Zhang et al. 2003; Morsony 
et al. 2007; Mizuta et al. 2011) show the jet drilling through 
the stellar envelope, pushing material towards the sides, 
forming a hot cocoon (e.g. Aloy et al. 2000). The surround- 
ing gas then acts to collimate the jet. Zhang et al. (2003) 
extracted angular profiles of the local rest mass density, total 
energy flux and outflow Lorentz factor at certain radii from 
the simulations. The resulting profiles can in many cases be 
approximated as constant within a characteristic jet opening 
angle, then decreasing as power laws towards the jet edge. In 
this work we adopt a similar parametrization of the angular 
profile of the bulk jet Lorentz factor, with a characteristic 
jet opening angle, power law index and normalization as free 
profile parameters. 

We develop a model for photon propagation in the con- 
text of a steady, optically thick, axisymmetric jet with angle 
dependent electron number density, photon number density 
and bulk Lorentz factor. We compute the observed spectrum 
taking into account contributions from the entire emitting 
volume as seen by an observer located at any viewing angle. 
As an example solution, we consider fireball dynamics as a 
way to relate the angle dependent parameters, in combina- 
tion with the assumed angular Lorentz factor profile. We 
developed a Monte Carlo simulation, unique in its ability to 
calculate photon propagation in a non-spherical explosion. 



We use this simulation to analyze the importance of angular 
bulk photon diffusion. 

We show that for a large region in the parameter space, 
the low energy photon index is approximately a w — 1 
(where dN/dE oc E a ). This is similar to the average value 
observed in GRBs (Kaneko et al. 2006; Nava et al. 2011; 
Goldstein et al. 2012). Furthermore, we present analytical 
expressions for the important energies and photon indices 
of the observed spectrum as functions of the free model pa- 
rameters. We show that photospheric emission by itself can 
account for observations of GRB spectra below the peak 
energy without the need for energy dissipation below the 
photosphere or additional radiative processes, provided that 
the characteristic jet opening angle is not much larger than 
~ few/T or if the outflow is viewed off-axis. Although we 
focus on GRB parameters, the results are general and can 
readily be applied to any optically thick, relativistic outflow 
with lateral outflow properties such as active galactic nuclei. 

This paper is organized as follows. In §2 we develop the 
analytical expression for the observed spectrum. We quali- 
tatively explain the spectral features in §3 in terms of con- 
tributions from different jet regions. The Monte Carlo code 
is explained in §4, and in §5 we present both simulated and 
numerically integrated spectra for outflow parameters rel- 
evant to GRBs. We discuss model sensitivity and photon 
time delays in §6 and summarize our results in §7. 



2 MODEL: PHOTON SCATTERING IN A 
STEADY JET 

Motivated by the results of Zhang et al. (2003) , we choose 
to define the shape of the angular Lorentz factor profile with 
an equation of the form 

(r — r \2 _ (rp — r m j n ) . . 

( minj ~ (e/e^ + i ' { ) 

where Fo is the maximum value of the Lorentz factor, 8 S is 
the characteristic jet opening angle, p determines the gradi- 
ent of the profile and F m i n = 1.2 is the lowest value of the 
Lorentz factor, differing from unity for numerical reasons. 1 
Three angular regions can be identified; the inner jet, 
the outer jet and the envelope. While this separation is arti- 
ficial and a consequence of our chosen Lorentz factor profile, 
it is useful for understanding the resulting observed spectra. 
The inner jet is characterized by an approximately constant 
Lorentz factor, F = Fq- In the outer jet the Lorentz factor 
decreases approximately as a power law of the angle with 
index p, and in the envelope the Lorentz factor has an ap- 
proximately constant value of F w 1. A characteristic an- 
gle separating the inner and outer jet regions is 0j, where 
T(#j) « Fq/VZ- Similarly, 9 C « #jr ( 1 / p is a characteristic an- 



1 As shown in §3.3, the exact value of T m i n only affects the very 
low energy spectrum, many orders of magnitude below the ob- 
served peak energy. 
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Figure 1. An example Lorentz factor profile (Eq. 1). In the inner 
jet (6 < 0j) the Lorentz factor is approximately constant, T = To 
while in the outer jet (9j < 6 < 6 e ) the Lorentz factor decrease is 
approximately that of a power law with index —p. In the envelope 
(9 > 9 C ), the Lorentz factor is approximately constant with a 
value of r m i n = 1.2. The dashed vertical line indicates 6j, while 

1 It) 

the dash-dotted line shows the location of 9 C ss 6jT ' . For this 
figure T = 100, 9j = 0.01 and p = 2. 

gle separating the outer jet and the envelope regions. 2 An 
example Lorentz factor profile is shown in Figure 1. 

Physical properties of the jet are expressed in spheri- 
cal coordinates (r,8,(f>), with the polar axis aligned to the 
jet axis of symmetry (see Figure 2). Some properties, such 
as the Doppler boost, depend on the angle to the LOS. 
These are most easily expressed in spherical coordinates 
(r, #los, 4>los), with the polar axis aligned to the LOS. 
When viewing the jet head on, the two sets of coordinates 
coincide, and the radiative contributions from the jet are in- 
dependent of azimuthal angle due to symmetry around the 
LOS. Our jet profile contains four parameters (including the 
viewing angle), three more than the simplified scenario of 
spherical symmetry. 

2.1 The photosphere and the decoupling radius 

The photospheric radius, R p h, is defined as the radius which 
fulfills the following condition: The optical depth for a pho- 
ton that propagates from that radius towards the observer 
is equal to unity. For the parameter values considered in this 
work, the comoving temperature at all relevant radii is much 
lower than the electron rest mass energy, and so we use the 
Thomson cross section in describing the scattering process. 
Using the variable z = rcos#LOS as the spatial coordinate 
in the direction pointing towards the observer, the optical 
depth from z = z m \ u to z — > oo is r ray = n' e <TT'D~ 1 dz, 
where T> = [F(l — /3 cos ^los)] -1 is the Doppler boost, n' e is 
the comoving electron number density and <tt is the Thom- 
son cross section (Abramowicz et al. 1991; Pe'er 2008). The 

2 In the calculations below we define 9 C to be the angle where 
r(9 e ) = \/2r m i n , which more accurately describes the character- 
istic angle separating the outer jet and the envelope regions. This 

gives e a = 6j[(V2 + ljro/r^n] 1 /?. 




</>LOS 



Figure 2. Schematic view of the jet. The solid lines represent the 
characteristic jet opening angle, 6j, as measured from the jet axis 
of symmetry, indicated by the dash-dotted line. The dotted line 
shows 9 C , the characteristic angle separating the outer jet and the 
envelope. The dashed line is the LOS, at an angle 9 V to the jet 
axis. Note that in general, 0LOS ^ 9 + 9 v . 

photospheric radius is a function of angle to the LOS, as fol- 
lows from the definition of r ray . For a spherically symmetric 
outflow, due to symmetry, the shape of the photosphere is 
independent of viewing angle. In this work we consider jet- 
ted outflows, for which this is no longer true. 

Here, in addition, we define the decoupling radius, -Rdcp, 
as the radius from which the optical depth for a photon 
that propagates in the radial direction is equal to unity. 
The optical depth from r = r m i n to r — >• co is T r di 
| r(l — f3)n' e <7Tdr. The importance of -Rdcp is that this 
is the characteristic radius where photons and electrons fall 
out of thermal equillibrium. If the outflow has axisymmetric 
comoving electron density or bulk Lorentz factor, then Rdcp 
is a function of angle to the axis of symmetry. For photons 
propagating along the LOS (cos#los = 1), dz = dr and 

= Rdcp- 



2.2 Formation of the observed spectrum 

In order to calculate the observed spectrum we consider the 
following: while the photosphere is defined as R p h = *"(r ray = 
1), in reality photons can scatter anywhere where electrons 
exist. Since we consider an expanding plasma in a steady 
state, electrons are assumed to occupy the entire volume. 3 In 
order to calculate the observed spectrum one has to integrate 
the emissivity over the entire volume, where the emissivity 
in volume element dV is proportional to the probability of 
photons to make their last scattering in dV . 

We assume that the local, comoving photon energy dis- 
tribution is well described by a Planck spectrum with tem- 
perature equal to the local, comoving electron temperature. 4 
As we consider the photon-electron cross section to be en- 
ergy independent, it follows that the optical depth is also en- 
ergy independent. In this section we calculate the observed 
spectrum assuming that it depends only on the last scat- 
tering positions of the observed photons. We justify this as- 



3 Since the probability of scattering atr> i? p h is negligible, this 
is a good assumption at late times. 

4 While this is true only for r < -Rdcpi as shown in §3.2 this is 
a good assumption for most jet profiles considered in this work. 
We do not use this assumption in §4. 



Rph 
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sumption by simulations (§4) which consider full treatment 
of photon propagation below the photosphere. 

Each time a photon in the outflow scatters into a di- 
rection pointing towards the observer, it has the probability 
cxp(— r ray ) to escape the outflow. The rate of photons es- 
caping the outflow from volume dV into solid angle dfl v is 



dNcsc = exp (— T ray ) dQ. v dV 



(2) 



where dh sc /d£l v is the number of scatterings within time 
dt in the volume element dV for which the outgoing pho- 
ton direction points within solid angle dfi v towards the ob- 
server. Generally, dh sc /dQ, v and T ray are functions of the 
viewing angle. Due to relativistic beaming, it is convenient 
to express dh sc /dQ v in terms of comoving frame properties. 
From here on, primed quantities are evaluated in the local 
frame comoving with the outflow. The number of scatter- 
ings within the volume element per unit time is Lorentz 
invariant, since dV' = FdV, dt' = T~ 1 dt and the total num- 
ber of scatterings only involves counting. The solid angle 
transforms as dfl v — V~ 2 dQ,' v (Rybicki & Lightman 1979), 
hence we write dh BC /dQ, v = V 2 dh' sc /dQ' v . The volume ele- 
ment is dV — r 2 dQ,dr, where r is the radial coordinate and 
dQ. is a solid angle element normal to the radial direction. 
Since the observed photon flux at luminosity distance dh 
is (l/dl)dN CBC /dQ. v , we write the expression for the energy 
spectrum, 



dflt 



exp(-r ray ) < E 



'dE 



dttdr, (3) 



where dP/dE is the probability density distribution describ- 
ing the normalized, local photon spectrum within volume 
element dV and the integration is over the entire volume. 

We now seek the probability for a photon within vol- 
ume element dV to have a lab frame energy between E and 
E + dE. Under the assumption that photons at r <C -Rdcp 
have a comoving Planck distribution with temperature T' , 
the probability of a photon to have energy between E' and 
E' + dE' is dP oc (B,y(T')/E')dE', where B„<(T') is the 
Planck spectrum with temperature T' . Normalization gives 
/ °° E' 2 /[exp(E'/kT') - \]dE' = (fcT') 3 C(3) « 2.40(fcT') 3 
where £ is the Riemann zeta function. The observed pho- 
ton energy and plasma temperature are E = DE' and 
T ob = VT , respectively, and therefore 



dP 
dE 



1 E 2 

2.40(fcT° b ) 3 exp(£/fcT ob ) 



1' 



(4) 



where k is the Boltzmann 

nob 



constant. From Eq. 4, 
E(dP/dE) oc B„(T ob ) within volume element dV. As the 
integration is over the entire volume, it follows that the ob- 
served spectrum is a superposition of blackbodies. 



2.3 Angle dependent jet properties 

We summarize the jet interaction with the surrounding envi- 
ronment as an angle dependent mass outflow rate per solid 
angle, dM/dSl = dM(6)/dfl, an angle dependent photon 
emission rate per solid angle, dN^/dtl = gL/V 7 (#)/gK7 and 
an angle dependent Lorentz factor, Y = Y(6), all symmetric 



around the jet axis. We further assume the Lorentz factor 
and mass outflow rates to be independent of radius. We omit 
writing angular dependences in the equations below for clar- 
ity. Assuming isotropic scattering in the comoving frame, 



dh' BC 

dnl 



ot cn e n 7 

4-7T 



(5) 



where c is the speed of light. The comoving electron and 
photon number densities are given by 



dM 



r 2 m p cl3r dfi 



and 



1 dN~ 



(6) 



(7) 



7 r 2 cT dQ 

respectively, where m p is the proton mass. Using Eq. 6, the 
decoupling radius can be written as 



Rdcp 



ctt dM 



(1 + P)PT 2 m p c dfi 



(8) 



where use was made of the assumed radial independence of 
T and dM/dQ. 

The photospheric radius is calculated as follows. Using 
the fact that rsin# L os is constant along a ray parallel to 
the LOS, dz = — r 2 (r sin^LOs) ^d^LOS- Using this change 
of variables, we write r ray = (<jt /m p cr sin#LOs) / L ° S [(^ ~ 
/3cosfeos)//3](dM/df2)d#LOS, from which we obtain 



R 



ph 



m p c sin 6> LO s 



e LOS 



/3 cos <9 L os dM 

8 In ' 



7LOS- 



(9) 



We thus find that for non-zero viewing angles, the pho- 
tospheric radius is a function of two variables (R p h = 
R P h(0LOS, 0los))- Since r ray oc r~\ we write r ray = R p h/r. 
Using Eqs. 3, 5, 6, 7 and 8 we obtain the observed spectrum, 



F° E h (8 v ) 



And 2 
Rdc 



jj + ^ 



dn 



■ exp 



-Rph^ 



'dE 



dttdr. 



(10) 



Eq. 10 (with use of Eqs. 4, 8 and 9) provides the 
energy flux for a given profile. It can be solved once 
T(6>), dM(6)/dQ, dN 1 {e)/dQ, and 8 V are known. In solv- 
ing Eq. 10 we consider non-dissipative GRB fireball dy- 
namics (Meszaros 2006) as a way to relate the angular jet 
properties. A similar treatment can be used in analyzing 
other outflow types. For simplicity, we assume an angle in- 
dependent luminosity, 5 dL/dQ, = L/4tt where L is the to- 
tal outflow luminosity. The photon emission rate from the 

5 As shown below, the part of the spectrum expected to be ob- 
served in a prompt GRB is formed by photons making their last 
scattering at approximately 9 < 5/ro- For model J A in Zhang 
et al. (2003), close to the largest radii of the simulation the 
Lorentz factor is To fa 140 up to 8 fa 0.02, while dL/dQ fa const 
up to 8 fa 0.05 fa 7/Fq (as seen in the top panels of Figures 8 and 
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base of the outflow, r — ro, is therefore also independent 
of angle, dN^/dQ = N-jAn and iV 7 = L/2.7kT where 
To = (L/A-Kroac) 1 / 4 is the temperature at the base of the 
outflow and a is the radiation constant. Any baryons and 
associated electrons present are rapidly accelerated by the 
intense photon field. As the outflow Lorentz factor grows, 
the angular separation required for causal contact between 
different parts of the outflow decreases. Therefore, we now 
make the assumption that a fluid element moving in a direc- 
tion within solid angle dfi evolves separately from the rest of 
the outflow, as if expanding spherically. The saturated value 
of the Lorentz factor is equal to the dimensionless entropy, 
77 = dL/dMc 2 , and so 

r= . (ft) 

4nc 2 dM/dQ, 

We choose to define T = T(6) as given by Eq. 1, and 
so dM{ff)/<Kl follows from Eq. ff. Due to adiabatic energy 
losses the comoving temperature above saturation decreases 
as T' = To(ro /r a )(r s /r) 2 / 3 , where r a = r/ro is the saturation 
radius. 



3 CONTRIBUTIONS FROM DIFFERENT JET 
REGIONS TO THE OBSERVED SPECTRUM 

The observed spectrum is obtained by solving Eq. fO, 
namely integrating the emissivity over the entire volume. 
While analytical integration is difficult for the different jet 
profiles, numerical integration is straight forward. In §3 - 3.3 
we assume the outflow to be viewed head-on. Therefore the 
radiative contributions from the different parts of the jet are 
symmetric around the LOS and Eq. 10 is readily integrated 
over azimuthal angle. Non-zero viewing angles are discussed 
in §3.4 and further explored numerically in §5. 

In Figure 3 we present the numerically integrated spec- 
trum where we have used the Lorentz factor profile defined 
in Eq. 1. For this figure we use the parameters T = 100, 
6j — 0.03 and p = 4, as well as 9 V = 0. The spectrum has a 
wavy shape, with a slope of roughly Fe oc E° over approx- 
imately 6 orders of magnitude in energy, which does not 
resemble the Planck spectrum. We identify several charac- 
teristic photon energies in the spectrum (Speak, E s , E-p, E c 
and .Eenv)- Analytical expressions for these energies are given 
below. 

The spectrum can be understood as a combination of 
components originating from different jet regions. Spectra 
integrated separately between the angular limits for the in- 
ner jet, outer jet and envelope regions are shown in Figure 
3 as dashed, dash-dotted and dotted lines respectively. The 
advantage of considering the emission from different jet re- 
gions separately is that it allows for analytical expressions 
of key spectral features. 

While the integration is over the entire volume, in or- 
der to explain the spectrum one can define characteristic 
photon energies as a function of angle to the jet axis. This 
is done in the following way. Consider only photons that 
make their last scattering between angles 9 and 9 + d0. The 

9 in Zhang et al. (2003)). Therefore, dL/dil = const is a good 
approximation for similar outflows. 
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Figure 3. The observed spectrum from a relativistic, optically 
thick outflow obtained by numerical integration of Eq. 10. Zero 
viewing angle is assumed. Separate integration of the contribu- 
tions from the inner jet, outer jet and envelope is shown with 
dashed, dash-dotted and dotted lines respectively. Photon ener- 
gies associated with characteristic outflow angles arc indicated by 
small vertical lines. For this figure To = 100, 8j = 0.03 and p = 4. 
A total outflow luminosity of L = 10 52 ergs -1 and base outflow 
radius ro = 10 s are used. 

distribution of last scattering radii spans a large range of 
radii (see Figure 4, right panel). At r <C R p h the flux is 
attenuated by the large optical depth, while at r > R p h 
the probability of scattering is low due to the optical depth 
being much smaller than unity. Therefore, the character- 
istic radius of the radial distribution is R p h(9). The dis- 
tribution of observed energies of photons that make their 
last scattering between angles 9 and 9 + d9 peaks at E w 
2.7/cT ob = 2.7V(9,r)kT'(9,R ph (9)). Therefore, one may 
consider E(9) = 2.7kT ob (9, R ph (9)) to be the characteris- 
tic energy of photons making their last scattering between 
angles 9 and 9 + d9. Typically, the Doppler boost decreases 
with angle while the photospheric radius increases with an- 
gle. 6 Therefore, the characteristic photon energy is expected 
to decrease with angle. Below we use this approximation to 
obtain the characteristic photon energies shown in Figure 3. 
For 9 V > 0, this discussion has to be generalized to include 
azimuthal angle. 

3.1 Inner jet component 

Emission from the inner jet region (0 ^ 9 ^ 9- s ) is shown 
with a dashed line in Figure 3. The most energetic part of the 
spectrum originates from this region due to the high Doppler 
boost and low photospheric radius (see Figure 4). Within 
this region the Lorentz factor is approximately constant. 
Therefore the observed spectrum from the inner jet region 
is similar to that of an optically thick, spherically symmet- 
ric wind with Y w Fo when only considering the emissivity 
within < 9 sC 9 S . The Doppler boost, V « 2Y/(\ + F 2 9 2 ), 

6 The exact angular behaviour of the Doppler boost and the pho- 
tospheric radius depends on the angular dependence of the out- 
flow parameters. 
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Figure 4. The left panel shows the photosphcric radius (-R p h, 
black) and the photon decoupling radius (-Rdcpi gray) as func- 
tions of angle. The dashed lines are the power law approximations 
used in §3.2. The right panel shows the same photospheric radius 
(black), and Monte Carlo simulations of the photon last scat- 
tering positions before reaching the observer (gray) for 8 V = 
(0.000 sg 8 V s£ 0.006). The Lorentz factor profile r = 100, 
8j = 0.03, p = 4 was used along with 8 V = 0. A total outflow 



luminosity of L 
radii. 



10 52 erg s 1 was assumed in calculating the 



and the photospheric radius for a spherically symmetric 
wind, R ph oc (1 + F 2 6 2 /3) (Pe'er 2008), are approximately 
constant for angles up to 9 ~ 1/To, and so the character- 
istic energy of photons making their last scattering within 
< 9 < l/r is approximately constant and is equal to 

-Epcak . 

The observed temperature of photons originating from 
angles larger than 1/ro decreases with angle due to the in- 
creasing photospheric radius and decreasing Doppler boost. 
For l/r < 9 < <9j, R ph /R dcp oc (l + r§6> 2 /3) (see Figure 4). 7 
Therefore, the approximation of thermal equilibrium close 
to the photosphere becomes less accurate with increasing 
angle (within the inner jet region). Since the characteris- 
tic photon energy is a monotonically decreasing function 
of angle for a spherically symmetric outflow, the approxi- 
mation of thermal equilibrium becomes less accurate with 
decreasing observed photon energy. By comparing numeri- 
cally integrated spectra to the full Monte Carlo simulations 
we find that for E > 10 _2 -E pe ak the assumption of thermal 
equilibrium can be used. 

The superposition of comoving spectra causes softening 
of the Rayleigh- Jeans part of the observed spectrum from 
the inner jet. Through simulations we find that the spec- 
tral index at E = 10 _2 -E pca k from a spherically symmetric 
outflow is approximately 0.6 units softer than the Planck 
spectrum. This agrees with the index found by Beloborodov 
(2010) who considered radiative transfer in a spherically 
symmetric, collisionally heated outflow. 

The peak energy of the inner jet component is ob- 
tained by considering the Doppler boosted comoving tem- 



perature at the photospheric radius along the LOS. It is 
the same as that from a spherically symmetric fireball with 
T = r , for which V « 2r , T' = (T /r )(r s /r) 2/3 and 
Rph = Lot /87rm p c 3 ro under the condition i? ph > r s (e.g. 
Meszaros 2006); 



J pcak 



5.4 x k 



Anr 2 , ac 

-5/12 



1/4 



/ s^mp^ro 



,3^4 \ 2/3 



ro 



= 540 



( L Y"' (ro \ 1/6 f r \ 

\w 52 J vio 8 / V 300 / 



8/3 



keV. (12) 



Here and below we use cgs units. All characteristic energies 
below are presented in units of -E pca k- At angles smaller than 
9j the Lorentz factor is approximately constant. Therefore, 
the observed spectrum is approximately that of a spherical 
wind for E > Ej = E(9j). 8 Under the assumption of thermal 
equilibrium at R p h(9j) we obtain the lower limit 



£j/£poak > 1 + 



r*ef 



i + 



-2/3 



(13) 



where the approximations sin( 



9j and j3j 



1/2T 2 were used. This expression simplifies to Ej/E pcai{ > 
4.2(r 6»j)" 10/3 for 8j > l/r . We thus conclude that observa- 
tions of spectra down to two orders of magnitude below the 
peak will reveal deviations from the spherically symmetric 
scenario if 9 } < 5/Fq- 



3.2 Outer jet component 

The spectrum from the outer jet region (0j < 9 < 9 C ) is 
shown with a dash-dotted line in Figure 3. It forms an ap- 
proximate power law between two limiting energies, E e and 
Et>- We derive analytical expressions for the power law in- 
dex and the upper and lower energy limits below. No simple 
analytical expression exists for the spectrum within the re- 
maining energy range Et> < E < E s . However, this energy 
range is typically small for the parameter space considered 
in this work. 

To simplify the calculation, we approximate the pho- 
ton energy distribution as a delta function centered on 
the average observed photon energy within volume ele- 
ment dV: dP/dE = 6(E - 2.7fcT ob ). This approach is 
validated by the numerical integration of the spectrum. 
Furthermore, we assume the outflow to be in thermal 
equilibrium at r = R p h(9). This assumption is justified 
in the outer jet region for p > 1 below. We therefore 
limit the analytical calculation to outflows with p > 1. 
The characteristic photon energy within volume element 
dV is 2.7fcT ob = (2.7DfcT /r)(r s /r) 2 / 3 . The delta func- 
tion variable is changed to r using the relation 5(E) = 
\dEfdr\- 1 5(r), 



E 



dP 
dE 



ES E- 2.7kT° 



-rS 



3/2N 



(14) 



7 This follows from R p ^ = R dcp at 8 = 0LOS = 0, while i?<jc P = 
const and i? p h oc (1+F 2 8 2 /3) for a spherically symmetric outflow, 
which is a good approximation for 9 < 8 S . 



8 In practice, the spectrum starts deviating at slightly higher 
energies since the Lorentz factor is not exactly constant within 
the inner jet region (sec Figure 3). 



© 2012 RAS, MNRAS 000, 1-15 



A theory of photospheric emission from collimated outflows 7 



where e = E/2.7kTo is the observed photon energy in units 
of the average photon energy at the base of the outflow. 
Using Eq. 14 one carries the radial integral in Eq. 10 which 
becomes 



N-, 3 
4ndl 2 



J IYo 1 V J 



exp 



sin0d0. 



(15) 



where 1 + /3 « 2, r s = Fro and dN 7 /dfl = N 7 /4tt were used. 

In order to solve the angular integral in Eq. 15, V, T>, 
Rdcp and i? p h are approximated as power laws of 9. In the 
outer jet region the Lorentz factor is well approximated as 



(16) 



For values of p > 1, the Doppler boost, V « 2r/(l + F 2 6> 2 ) 
has a break at 0t> = 9 i {F 9 i ) 1/(l '- 1) (see Figure 5). 9 For 
angles 9t> < 9 < 9 C the desired parameters can be expressed 
as power laws in to a good approximation, and so we solve 
for the power law spectrum between these two integration 
limits. The characteristic photon energies associated with 
6-d and 9 C gives the upper and lower energy limits Et> and 
E e . For 8t> < < 9 C the Doppler boost is approximately 



V = 2T Q 



(17) 



Using the approximate Lorentz factor profile in Eq. 16 the 
decoupling radius (Eq. 8) is readily obtained, 



^Idcp — p 3 
1 



3p 



(18) 



where _R* = Lot /8ivm p c 3 and /3 « 1 were used. Us- 
ing the small angle approximation sin 9 « 9, the photo- 
spheric radius (Eq. 9) at angle 9 is approximately R p h — 
(2R*/0) / * r~ 2 £> _1 cfl9. The integral may be split at angle 
9t>, representing contributions to the optical depth from 
electrons at larger and smaller angles, respectively. The con- 
tribution from within 9t> is subdominant for p > 1 and 
9 ^> 9t>, and so it is sufficient to consider 9t> as the lower 
angular limit. Using Eqs. 16 and 17 the angular integration 
for the photospheric radius is readily solved, 



R P h — 



R, 



(3 P + i)r 3 



3p 



(19) 



which is a factor 3p + 1 lower than Rd cp for all angles 9t> < 
9 < 9 C . We therefore conclude that the outflow is in thermal 
equillibrium at the photosphere within this angular range. 
The approximated power laws are shown in Figures 4 (left 
panel) and 5. 

In order to find the photon energy associated with 
9-d, E v » 2.7V(8 T ,)(kT'/r(9 T ,))(r s (9 T ,)/R ph (9 T ,)) 2 ' 3 , one 

9 The break is located at 6t>T(0d) = 1, and Y is given by Eq. 16. 
For p = 1, Y 2 9 2 = const and there is no such break. 



first calculates R p h(9i>). 10 Using this approximation and 
V(6t>) = F(6t>), we calculate the photon energy associated 
with 9t>, 



(p + 3) 2/3 
Et)/E pc ^ = (ro^j) 



8p/3(p-l) 



(20) 



which is valid for 8j > 1/Tq. As can be seen from the def- 
inition of 9t> (below Eq. 16), 9t> ~ #j when 9j w 1/ro. 
Therefore Et> « -Ej ~ Speak when 9j w 1/Fo. Thus, for nar- 
row jets (9j « 1/Fo) the upper limiting energy to the power 
law is approximately the same as the peak from the inner 
jet region. For 9j > 1/Fo there is a bump-like spectral fea- 
ture related to Ed (see Figure 3). If this feature is observed, 
constrains can be put on the jet properties through Eq. 20. 

Inserting Eqs. 16, 17, 18 and 19 into Eq. 15 gives the 
observed flux within the energy range E e <C E <C Et>, 



^ob iV 7 3 4 R, /e\3/2 

f K 



Aivd 2 2 (F 9 p ) 2 r 



1 1 3p+l (r flf) 4 r \2) 
N 1 3 (3p + 1)K? +1 ) ^ a (R t \ 5( 1_ p) 



(Fo0fP 



4-7rd 2 2 p 1 " \r 

(0 4< P> y x 3 (p _1 ) exp(— a;)da; 



(21) 



where x = (3p + l)" 1 ^? )-*(R./r )(e/2) 3 ' 2 e ip , x min = 
(3p + l)- 1 (R t /r () )(e/2) 3 ' 2 (r 9f) 4 ^- 1 and z max = (3p + 
l)- 1 ( J R»/r )(e/2) 3 / 2 ((l + ^)/r min ) 4 . 

The integral in Eq. 21 is approximately constant for a 
large range of energies. The integrand decreases exponen- 
tially for large x. When a; max S> 1, the upper limit is well 
approximated as infinity. This condition translates to 



Ec _ I ,„ , t\ r ( Tmin 

e>> 2TMo= 2 \ (3p + 1) MT77I 

In units of the peak energy, one obtains 



2/3 



E c j Speak 



(3p + i) 2 / 3 /r min 



(V2 + 1) 8 / 3 V To 



8/3 



(22) 



(23) 



Setting the upper limit to infinity, the integral in Eq. 21 
becomes 



oc 

J x^ip' 1 ) exp (-a;) da; = r inc Q + 1^ , x min | 



(24) 



10 At 9d, the electrons at angles 9 < 6-q have to be taken into 
account when calculating the optical depth. By integrating the 
expression for the photospheric radius (Eq. 9) over angles 9 < 9-p, 
one obtains R ph a (R»/r 3 )(r 6»j)- 1 /(P- 1 ) {1 + (r o 0j) 2 [1/3 + (p + 
3)- 1 (r e i )(f+ 3 )/(f- 1 )(l - (T 9 i )-(P+ s y(P- 1 ))]}. This simplifies 
to Rp b ^ (R»/r 3 )(p + 3)- 1 (r 6 i ) 3p K p - 1 '> for 6>j > l/r . Note 
that this expression is different from Eq. 19, which is derived for 
9 > 6-d. 
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Figure 5. The Dopplcr boost (solid black) and the Lorcntz factor 
(dashed black) as functions of angle to the jet axis. 9 V = is 
assumed. The respective approximations used in §3.2 are shown 
as gray lines. The Doppler boost has a bump at angle 9 « 6x>- 
The characteristic photon energy associated with 9-p is shown in 
Figure 3. A Lorcntz factor profile (Eq. 1) with parameter values 
T = 400, 0j = 3/r , p = 4 was used. 

where Fi nc {(l/p+l)/2, a; m i n } is the incomplete gamma func- 
tion, which for a; m i n <§C 1 and p > 1 is approximately con- 
stant. 11 The condition x m in -C 1 translates to 

6 « 2 {(3p + 1) g (r efy^ } 2/3 « |g. (25) 

We thus find that the integral in Eq. 21 is approximately 
constant within the energy range E e <C E <C Ex>. The ob- 
served flux in this energy range is a power law approximately 
satisfying 

F£ oi E^-p) (26) 

for p > 1. Expressed as the photon index within this energy 
range, Eq. 26 implies a « — (1/4)(1 + 3/p). Redoing the 
above calculation for p = 1 results in the same power law 
index, with the modification that the upper limiting energy 
is set by Ej, and so Eq. 26 can be considered valid for p ^ 1. 
We therefore conclude that for a Lorentz factor profile with 
p = 1, a « — 1 while for p = 2, a ~ —0.5. For an angular 
profile with p = 4 and #j = 3/To (as shown in Figure 3), 
■Ed /Speak = 2.5 x 1CT 2 and E e /E pcak = 4.6 x 1CT 5 . For jets 
with large opening angles (9 S > 10/r ( )), the range of energies 
between Et> and E c is small and so the outer jet spectrum 
is not well approximated as a power law. 

3.3 Envelope component 

The spectrum from the envelope component (9 < 9 C ) is 
shown with a dotted line in Figure 3. It is characterized 

11 For •^min — the incomplete gamma function equals the 
complete gamma function, T{(l/p + l)/2}. For 1 < p < 10 2 , 
1 <r{(l/p + l)/2} < 1.76. 



by a peak at E — E anv , where E cnv is given by the observed 
temperature of the envelope component. 

The Lorentz factor of the envelope is approximately 
constant, close to unity. Due to the low bulk speed, the 
photospheric radius of the envelope is very large. There- 
fore, the comoving photon temperature at transparency is 
very low. Furthermore, the observed envelope temperature 
is barely Doppler boosted. Therefore, its spectrum is ap- 
proximately that of a blackbody with observed temperature 
equal to the comoving temperature at the radius of trans- 
parency. For r m ] n = 1.2, the envelope speed is /3 m i n ss 0.5, 
while cos 9 C « 1 for most model parameters. Therefore, the 
Doppler boost is V w [r min (l - /3 m i n )] _1 = r min (l + /3 min ). 
The photospheric radius for angles larger than 9 e is approx- 
imately constant, equal to the decoupling radius since the 
electrons move at a bulk speed that is only mildly relativistic 
(see Figure 4). Thus we calculate the ratio of the character- 
istic observed envelope photon energy to the peak photon 
energy, 

Eenv/iW - ( 1± |^) J/3 ^ ("ff ) 8/ '- (27) 

For an angular profile with To = 100 and F m i n = 1.2, 

-Eenv/Epeak ~ 3.3 X 10 



3.4 Non-zero viewing angles 

Due to the complexity of the calculation in this scenario, we 
give here only a qualitative explanation. Full quantitative 
treatment is presented in §5. Increasing the viewing angle 
has two main consequences: decreasing E pcl> ^ (and therefore 
also the bolometric flux) and softening the spectrum. Since 
the characteristic angles for the inner jet are smaller than 
those for the outer jet, the inner jet spectral component is 
aflected more than the outer jet component by viewing angle 
variations. 

As long as 9 V <JC 9j the observed spectrum is well ap- 
proximated by the 9 V — solution. For 9 V > 9- s , i5 p0 ak de- 
creases due to the decreasing Doppler boost as well as the 
increasing photospheric radius. Photons from the outer jet 
region contribute to a general softening of the spectrum be- 
low the peak energy. For 9 V > 9,, the total observed flux 
decreases with increasing viewing angle, and at some view- 
ing angle (dependent on jet parameters as well as detector 
characteristics) the flux in non-detectable. 

For p — > oo the jet profile approaches a top-hat shape. 
For such jets the peak energy and bolometric flux decreases 
rapidly with increasing viewing angle (for 9 V > 9,) and the 
probability of observing the outflow at 9 V > 9j is low. How- 
ever, for outflows with a shallow Lorentz factor gradient 
(small values of p) , the bolometric flux and peak energy de- 
crease are slower. Therefore the outflow may be observed up 
to # v w few x 9j , depending on detector characteristics and 
the exact jet profile. Assume # Vimax to be the largest view- 
ing angle where the outflow is still detectable. The fraction 
of GRBs observed at V < 9j is Afij/An v , max = 9 2 J9l >rn ^. 
For #v,max > \/2 x 6>j , the majority of observed outflows are 
viewed at 9 V > 9 y Therefore, the consequences of observing 
the outflow at 9 V w 9 S should not be neglected. 
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4 NUMERICAL SIMULATION 

In order to validate the analytical model and to explore the 
importance of full photon propagation history, a numerical 
code has been developed. The code is based on an earlier 
code for photon propagation in a spherically symmetric, rel- 
ativistically expanding plasma (Pe'er & Waxman 2004; Pe'er 
et al. 2006; Pe'er 2008). This code is specially designed to 
allow for the bulk Lorentz factor, mass outflow rate and 
photon emission rate to be functions of angle with respect 
to the jet axis. This is required in the current context as the 
optical depth between two points in the outflow depends 
on the local Lorentz factor and electron density. As pho- 
tons propagate across the jet, photon energy changes due to 
propagation between regions of different bulk Lorentz fac- 
tor is inherently considered; therefore the comoving photon 
spectrum at the last scattering position is not artificially 
fixed to blackbody. No assumption of thermal equillibrium 
is made above the photon injection position, as opposed to 
the analytical model. The full Klein-Nishina cross section is 
used in the scattering process. 



4.1 Code description 

The code tracks photon propagation, starting from a ran- 
dom position in the outflow with T r( ji = 20, ensuring a low 
probability (exp(— 20) ~ 2.1 x 10~ 9 ) for a photon to escape 
the outflow without being scattered. The lab frame angu- 
lar coordinates of the initial position are randomly drawn 
from a uniform distribution. The initial photon propaga- 
tion direction is randomly chosen in a similar way. The 
initial comoving photon energy is drawn from a black- 
body distribution with temperature equal to the comov- 
ing outflow temperature at the photon position; T'(6,r) — 
T (r /r s (9))(r s (9)/r) 2 ^. 

A photon has a probability exp(— r) to propagate an 
optical depth r before scattering. Thus, an optical depth, 
At, representing the distance between scatterings in the di- 
rection of propagation, is drawn from a logarithmic distri- 
bution. The optical depth from the photon position to in- 
finity in the photon propagation direction, Too , is calculated 
and compared to the drawn optical depth. If At ^ Too the 
photon escapes the outflow, else the distance corresponding 
to Ar is calculated, taking into account the angular depen- 
dence of the Lorentz factor and electron number density. The 
new photon scattering position is then considered. The pho- 
ton four-vector is transformed to the local comoving outflow 
frame. The scattering electron Lorentz factor is drawn from 
a Maxwellian distribution with temperature T'(8,r) along 
with a random electron propagation direction. The photon 
four- vector is transformed to the electron rest frame and the 
photon undergoes a Compton scattering, changing its energy 
and direction. The photon four-vector is transformed back 
to the lab frame, and a new optical depth is drawn. The pro- 
cess is repeated until the photon has escaped the outflow. 

After simulating a sufficient number of photons (typi- 
cally 3 x 10 7 ) the escaping photons are binned in viewing 
angles. The bin width is dependent on how sensitive the 
spectrum from a given profile is for viewing angle variations. 
The bin width used for each spectrum is presented in the re- 
spective figure caption. A typical width is A8 V ~ 0.005. 




Obs. energy, E/m c c 



Figure 6. Simulated (colored) and numerically integrated (black) 
spectra for a narrow jet observed at different viewing angles. In 
this plot, Tq = 100, 8j = l/To and p = 1. Three different viewing 
angles are shown, 8 V = (0.000 8 V ^ 0.0045, red diamonds and 
solid black lines), 8 V = 8j (0.009 ^ 8 V ^ 0.011, green triangles and 
black dashed lines) and 8 V = 2 X 8j (0.019 ^ 8 V ^ 0.020, magenta 
squares and dash-dotted lines). After viewing angle binning, the 
red, green and magenta spectra contain 1002, 1758 and 1486 pho- 
tons, respectively. The photon index below -E pGa k is a = — 1 for 
all viewing angles, two units less than the Rayleigh-Jeans index. 
For this parameter space region the numerical integration gives 
an excellent fit to the simulated spectra. 



5 RESULTS FROM SIMULATION AND 
NUMERICAL INTEGRATION OF THE 
SPECTRUM 

The obtained spectra are presented in Figures 6 to 10. 
In each figure we present both the simulated spectra and 
the numerical integration of Eq. 10. We have considered a 
large set of the parameter space region: To = 100, 200, 400; 
p = 1, 2, 4; 6>jr = 1, 3, 10 and 8 v /9 i = 0, 0.5, 1, 1.5, 2. 
In Figures 6 to 10 we have assumed non-dissipative fireball 
dynamics with L — 10 52 ergs _1 and ro = 10 s cm. A lumi- 
nosity distance of = 4.85 x 10 28 cm (z — 2) was used for 
spectrum normalization. Each run contains 3 x 10 7 photons 
(unbinned) unless otherwise noted. 

For a top-hat jet (p — > oc), the most likely viewing angle 
is 9 v /0j w 2/3 since the jet is obscured for 9 V > By However, 
the outflows we consider here can be observed at 8 V > 9j due 
to the gradual decrease of the Lorentz factor. This increases 
the value of the most likely viewing angle. A value of 9 V « 9j 
may be close to the average value for jets with intermediate 
values of p (the exact value of the most likely viewing angle 
depends both on the outflow as well as detector properties). 
Therefore we choose to show 9 V /9 S — 0, 1, 2 in Figures 6 
and 8. 

The most important result of this paper is that we ob- 
tain (irj-1 (where a is the low energy photon index in the 
Band function, dN/dE oc E a below -E pea k) for a large region 
of the considered parameter space. This is demonstrated in 
Figures 6 to 10. We consider this a very important result as 
this is similar to the average low energy photon index ob- 
served in GRBs (e.g. Kaneko et al. 2006; Nava et al. 2011; 
Goldstein et al. 2012). 
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Figure 7. Simulated (colored) and numerically integrated (black) 
spectra for narrow jets with different Lorentz factor gradients. In 
this plot, T = 100, 0j = l/r and 8 V = 6 i (0.009 < 8 V < 0.011). 
Three different values of p are shown, p = 1 (red diamonds and 
solid black lines), p = 2 (green triangles and black dashed lines) 
and p = 4 (magenta squares and dash-dotted lines). After viewing 
angle binning, the red, green and magenta spectra contain 1758, 
1900 and 891 photons, respectively. The low energy photon index 
is close toaRJ —1 for all Lorentz factor gradients considered here. 
For p = 4, the high energy spectrum does not decay exponentially 
due to photon diffusion from high angles. See further discussion 
in the text. 

For a narrow jet (9j — 1/Tq, see Figure 6), the peaks of 
the emission from the inner and outer jet regions (Speak and 
Et> respectively, see discussion in §3) coincides, and so the 
spectrum below iJp C ak is a pure power law for ~ 6 decades in 
energy. The resulting spectral shape is independent of view- 
ing angle, however the peak energy decrease as the outflow 
is observed at large viewing angles. 

In Figure 7 we present spectra for narrow jets with 
p = 1, 2 and 4 observed at 8 V = 6j. As can be seen, the 
low energy photon index is not very sensitive to the Lorentz 
factor gradient. The softest spectrum is obtained for p = 1, 
as shown in Figures 6 and 7. For higher Lorentz factor gradi- 
ents there is a slight increase in the photon index, consistent 
with the analytical expression in Eq. 26. For 1 < p < 4, the 
photon index is —1 < a < —0.5. Values of p lower than 
unity has been considered. As p decreases below unity the 
photon index below the peak energy increases. For p — 
the outflow profile is spherically symmetric. Therefore the 
observed spectrum is that of a spherically symmetric wind 
(a ~ 0.4). 12 For the simulated spectra, the exponential cut- 
off expected at energies above the peak energy becomes less 
sharp for increasing values of p. We discuss this further be- 
low. 

The spectra from jets with large opening angles (#j m 
10/To) observed at 6 V <C 6*j appear as those from spherically 
symmetric winds. However, for viewing angles 8 V ~ 6j the 
observed photon index below the peak energy is lower. In 

12 For values of p < 1, the angle separating the outer jet region 
and the envelope (8 C 8jF^ p ) becomes larger than unity (for 
8j ^ To) and the outflow consists only of the inner and outer jet 
regions. 



Figure 8. Simulated (colored) and numerically integrated (black) 
spectra for a wide jet observed at different viewing angles. In this 
plot, To = 100, 0j = 10/r and p = 4. Three different viewing 
angles are shown, 8 V = (0.000 8 V ^ 0.006, red diamonds and 
solid black lines), 8 V = 8j (0.009 ^ 8 V s£ 0.011, green triangles 
and black dashed lines) and V = 2 X 8j (0.019 s£ 6» v s£ 0.020, 
magenta squares and dash-dotted lines). After viewing angle bin- 
ning, the red, green and magenta spectra contain 1577, 1826 and 
1112 photons, respectively. For 9 V 9j, the high energy spec- 
trum resembles that of a spherical wind. For 8 V = 8j, the photon 
index is a = —1, similar to the spectrum from a narrow jet shown 
in Figure 6. Due to the large Lorentz factor gradient (p = 4), for 
8 V = 2 X 8j the peak energy is very low. 



Figure 8 we present spectra from an outflow with the profile 
To = 100, 9j =0.1 and p = 4 observed at different viewing 
angles. For 9 V = 0j, a m — 1 just as for narrow jets. Due 
to the large Lorentz factor gradient, i? pea k decreases rapidly 
with increasing viewing angle. In Figure 8, -Ep ea k(#v = 2 x 
S})/ E P cnk(dv = 0) ~ 10~ 3 . As discussed above, depending 
on the jet properties and detector characteristics the most 
likely viewing angle may be close to the jet opening angle. 

In Figure 9 we present observed spectra from jets with 
different opening angles (#jTo = 1,3 and 10) viewed head on. 
Both 6jTo = 1,3 result in low energy slopes close to a = — 1 
independent of viewing angle, while wider jets viewed at 
8 V — results in a spherically symmetric spectral shape a 
few decades below -E p0 ak- 

Increasing the maximum Lorentz factor, Fo, shifts the 
spectral components from the inner and outer jet up in en- 
ergy while keeping their spectral shapes intact (Figure 10, 
also see Eqs. 12, 13, 20 and 23). The envelope component 
is unaffected (as expected, see Eq. 27). As long as To#j is 
constant, the spectral shapes of the inner and outer jet com- 
ponents are not affected by varying To. 

We reach the conclusion that the low energy spectral 
index is close to a ~ — 1 for narrow jets, with 6*j ^ few/To 
and moderate Lorentz factor gradients for all viewing angles. 
Similar photon indices are obtained from wider jets observed 
at 8 V w 9j . 
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Figure 9. Simulated (colored) and numerically integrated (black) 
spectra for jets with different jet opening angles. In this plot, 
T = 100, p = 2 and 8 V = (0.000 < 8 V sg 0.006 for the red 
and green spectra, 0.000 ^ 8 V ^ 0.020 for the magenta spectra). 
Three different values of 8j are shown, 9j = I/Fq (red diamonds 
and solid black lines), 9 S = 3/To (green triangles and black dashed 
lines) and 8j = 10/To (magenta squares and dash-dotted lines). 
After viewing angle binning, the red, green and magenta spec- 
tra contain 2163, 2267 and 4482 photons, respectively. The exact 
spectral shape depends on the value of p, however the average 
spectral index is approximately —1 for jet opening angles up to 
6j Ri few /Vq and values of p up to a few. 



be 
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Figure 10. Simulated (colored) and numerically integrated 
(black) spectra for narrow jets with different values of To but 
constant values of djT'o- In this plot, p = 1, 6>j = 1/To and 9 V = 
(0.000 < 9 V < 0.006 for red, 0.0000 ^ 8 V ^ 0.0032 for green and 
0.000 ig 6» v ^ 0.0020 for magenta). Three different values of T are 
shown, To = 100 (red diamonds and solid black lines), To = 200 
(green triangles and black dashed lines) and To = 400 (magenta 
squares and dash-dotted lines). After viewing angle binning, the 
red, green and magenta spectra contain 1934, 558 and 1114 pho- 
tons, respectively. The magenta spectrum requires small viewing 
angle bins due to the high To- Therefore, 1.3 X 10 8 photons were 
simulated for the magenta spectrum. For the other two spectra 
1 X 10 7 photons were simulated. The spectral components from 
the inner and outer jet region are shifted to higher energies in 
accordance with Eqs. 12 and 20. 
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Figure 11. Same as Figure 7, but 8 V = (0.000 s£ 8 V ^ 0.006). 
After viewing angle binning, the red, green and magenta spectra 
contain 1934, 2166 and 1178 photons, respectively. A power law 
above the peak energy is seen in the magenta spectrum. The 
power law is further discussed in the text. 

5.1 Asymmetric photon diffusion and 
Comptonization 

Preferential photon diffusion towards the jet center is ex- 
pected to occur in outflows where the bulk Lorentz factor de- 
creases from the jet axis. This can be understood in terms of 
average photon scattering angles. The average photon scat- 
tering angle with respect to the radial direction is ~ 1/T 
in the lab frame. Therefore a photon is more likely to scat- 
ter from a region of low Lorentz factor to a region of high 
Lorentz factor than the other way around. The importance 
of photon diffusion is dependent on the Lorentz factor gra- 
dient, and so the spectra from jet profiles with large values 
of p are expected to differ from the numerically integrated 
spectra where we assume dN-y/dQ. is r-independent. This is 
indeed observed in the simulations. The effects of angular 
photon diffusion are, however, sub-dominant as compared 
to the observed temperature variations across the jet when 
considering the formation of the observed spectrum from jets 
with small Lorentz factor gradients (approximately p < 4). 

Comptonization of photons propagating to regions of 
higher Lorentz factor may cause photons to fall out of ther- 
mal equillibrium with the local electrons. This has been in- 
vestigated by comparing the simulated lab frame photon 
energy to the lab frame temperature of the electrons at the 
last scattering position for all photons in a given spectrum. 
The observed photon energies are well described by the local 
observed temperature for most of the considered parameter 
space. However, a deviation is observed in Figures 7 and 11. 
The expected high energy cut-off is significantly smoother 
for the spectrum with p = 4 observed at 8 V = 6j. The same 
jet observed at V = forms a high energy power law with 
photon index /3 ~ — 2. The photons forming the power law 
start out with thermal energies in the outer jet region where 
the Lorentz factor gradient is large. Repeated scatterings 
towards the jet center increases the photon energy until the 
photon reaches the inner jet region where it finally escapes. 
The final photon temperature depends on the initial ther- 
mal photon energy, the number of scatterings before escape 
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and the average energy gained in each scattering. The last 
two terms are dependent on the Lorentz factor gradient. 
The Comptonization of photons propagating in media with 
a Lorentz factor gradient will be further explored in future 
works. In general, the analytical expressions for character- 
istic energies and spectral shapes presented in §3 agree well 
with the Monte Carlo simulations except for large Lorentz 
factor gradients. 



6 DISCUSSION 

6.1 Sensitivity to the chosen dynamics and 
Lorentz factor profile 

In this work we consider fireball dynamics for the outflow 
properties. However, other dynamics are possible: kinetic 
energy dissipation may be significant during the acceleration 
phase, or the outflow may be dominated by magnetic fields. 
The overall shape of the spectral components from the inner 
jet and evelope regions are not expected to be sensitive to the 
exact dynamics, however the value of -E poa k is dependent on 
the radial scalings of the outflow properties. The exact value 
of the power law slope in the outer jet component is expected 
to be modified for other dynamics, as well as Ex> (Eq. 20, the 
energy of the spectral bump below the peak energy in Figure 
3). However, the ratio Ej/E pe& k which gives an indication 
of the energy where the spectrum starts to deviate from 
the sperically symmetric spectrum is not sensitive to the 
dynamics. 

The shape of the angular profile is motivated by hy- 
drodynamical simulations. One may consider other shapes 
(e.g. Gaussian) for the angular parameters. The exact values 
of the characteristic energies and spectral shapes are natu- 
rally dependent on the chosen profile, however the general 
spectral features are not sensitive as long as the considered 
profile has similar characteristics. 

6.2 Relative photon time delays 

In this work we focus on the spectrum from steady jets. In 
much the same way as photons from a certain angle has a 
characteristic energy, they also have a characteristic time de- 
lay as compared to a photon from the center of the jet. One 
can therefore associate a given photon energy in the spec- 
trum with a certain time delay. Assuming the jet is viewed 
at zero viewing angle, the characteristic lab frame time delay 
as a function of angle is 



At(6) » ^iM[i-^(0)cos0]- [l-/?(0)] 



c 

RM9) 



[1-/3(61) cos 6»] 



(28) 



where 6 > 1 /To is assumed in the last step. The time delay 
for a given energy is generally larger for an outflow with a 
Lorentz factor that decreases from the jet axis, as compared 
to the sperically symmetric case. It increases with decreas- 
ing 6j and increasing p, but decreases rapidly with To. As 
an upper limit, we consider a jet profile with 9 } = 1/Yq. 
The time delay of photons with lower energy than the ob- 
served peak energy is calculated as follows. Consider photons 
that make their last scattering in the outer jet region. For 



6j = 1/To, 8t> = 0j and Eqs. 16, 17 and 19 are good ap- 
proximations for r(#), T>(9) and R p h(0) respectively. The 
characteristic photon energy at angle 6 is -E(6 l )/i5p 0a k ~ 

{[V(e)/rm[rs(e)/RM0)} 2/:i }/{2[r s (0)/RM0)] 2/:i }- Since 
At « R ph (l - /3cos6>)/c = R ph /VTc, we write 
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which is valid for energies much lower than the peak energy. 
Time-resolved analysis of spectra detected down to two or- 
ders of magnitude below the peak energy may resolve the 
superposition of blackbodies in the lowest observed energies 
if the time bins are smaller than At ob . The result in Eq. 29 
shows that for a narrow jet with otherwise typical GRB pa- 
rameters, the expected time delay is w 1 s for photons with 
energy E « 10~ 2 -E pca k. This is in sharp contrast to naive 
1-D calculations used in the past, which gives time delays of 
only few tens of milliseconds for the same parameter values. 
The spectral component associated with the envelope has 
a very long associated time delay, due to the large photo- 
spheric radius (see Figure 5). We therefore do not expect the 
envelope component to be observed in the prompt emission 
emitted from a transient source. 



6.3 Band-like GRB spectra 

Observed GRB spectra do not, in general, appear thermal. 
The spectra are usually well fitted by the Band function, al- 
though more complex spectra has been observed (e.g. Ryde 
et al. 2010; Ackermann et al. 2011; Guiriec et al. 2011). 
The average values of the low and high energy photon in- 
dices are a « — 1 and /3 « —2.5 respectively. This is in 
contrast to a blackbody spectrum in which a = 1 and the 
spectrum cuts off exponentially above the peak energy. To 
reconcile the photospheric model with observations, energy 
dissipation close to the photosphere has been considered as a 
way to modify the local comoving photon spectrum (Thomp- 
son 1994; Spruit et al. 2001; Rees & Meszaros 2005; Pe'er 
et al. 2005; Giannios & Spruit 2005; Pe'er et al. 2006; Laz- 
zati et al. 2009; Beloborodov 2010; Ryde et al. 2011; Zhang 
& Yan 2011; Vurm et al. 2011). In these scenarios the ob- 
served spectrum results from a combination of radiative pro- 
cesses. Comptonization of the thermal photon spectrum is 
the source of the high energy power law, and synchrotron 
photons contribute to the low energy spectrum (e.g. Vurm 
et al. 2011). In this work we show that three dimensional 
effects arising from the jet shape may add a significant con- 
tribution to understanding the soft low energy spectrum. 
Narrow jets, or wider jets observed at 6 V « 6 l j result in a low 
energy photon index a « — 1 under natural assumptions. 
Harder values of the low energy photon index are obtained 
for wide jets viewed at 9 V <ti S . 

In some GRBs observed by the Fermi Gamma-Ray 
Space Telescope (GRB 100724B, Guiriec et al. 2011 and 
GRB 110721 A, The Fermi LAT Collaboration et al. 2012), 
a spectral "bump" is significantly detected at energies be- 
low the main peak. This bump is commonly interpreted as 
photospheric emission (e.g. Ryde 2005) . In such a scenario, 
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the Band component originates from non-thermal emission 
processes in the optically thin part of the outflow outside 
the photosphere (e.g. Zhang et al. 20f2; Veres et al. 20f2). 
Within this interpretation the analysis performed in this pa- 
per is applicable to the bump instead of the Band component 
and may explain bumps wider than blackbody. 



7 SUMMARY AND CONCLUSIONS 

In this work we develop the theory of photospheric emission 
from relativistic jets with angle dependent outflow proper- 
ties. We consider a three parameter angular Lorentz factor 
profile (Eq. I) where the Lorentz factor is approximately 
constant, F = r within a characteristic jet opening angle 
9j and then decreasing approximately as V oc 9~ p towards 
the outer jet edge. The shape of the profile is motivated by 
the results of hydrodynamical simulations of jet propagation 
through the progenitor star (e.g. Aloy et al. 2000; Zhang 
et al. 2003; Morsony et al. 2007; Mizuta et al. 2011). In §2 
the expression for the observed spectrum is obtained analyt- 
ically by integrating the local emissivity over all radii and 
angles (Eq. 10). We derive approximate analytical expres- 
sions for the important spectral features in §3. Comparing 
jetted outflows to spherical outflows, we show that soften- 
ing of the spectrum below the peak energy is expected for 
an on-axis observer. This is a consequence of weaker radial 
beaming of photons at high latitudes due to the decreasing 
Lorentz factor. A Monte Carlo simulation was developed to 
investigate the importance of full photon propagation below 
the photosphere (§4). In §5 we present spectra obtained by 
numerical integration of Eq. 10 as well as the Monte Carlo 
simulations for different profile parameters and viewing an- 
gles (Figures 6 to 11). 

The most important result of this paper is that we ob- 
tain a photon index a £s — 1 below the peak energy for a 
large region of the considered parameter space. For narrow 
jets (dj < few /To) with Lorentz factor gradients 1 < p < 4 
observed at any viewing angle, we find —1 > a > —0.5. 
For jets with 9- 3 « 1/F and p ^ 1, a w -(1/4)(1 + 3/p) 
(see Eq. 26). Observing wider jets (9 S > 5 /To) at 9 V « 9j 
results in similar soft spectra with a « — 1. The spectrum 
from wide jets observed at 9 V <C 9 S appears similar to the 
spectrum from a spherical wind (Beloborodov 2010) for sev- 
eral decades below the peak energy. However, observing the 
outflow at viewing angles 8 V « is less likely than 9 V « 9 S . 
Additionally, increasing the viewing angle causes the ob- 
served peak energy to decrease. The decrease is slower for 
low Lorentz factor gradients. Therefore, such outflows may 
still be observed at 9 V « few x 9 S . 

Photon diffusion primarily towards regions of higher 
Lorentz factor is observed in the simulations. For outflows 
with large Lorentz factor gradient (p > 4) this bulk propaga- 
tion of photons has to be taken into account when computing 
the observed spectrum. Comptonization of the photons that 
make repeated scatterings towards regions of larger Lorentz 
factor can produce photons with energies significantly higher 
than the local temperature. Evidence for this can be seen in 
Figure 11 (p = 4), where an approximate high energy power 
law is formed above the peak energy. Photon propagation in 
plasma with a steep Lorentz factor gradient will be further 
explored in future works. 



The observed photospheric spectrum from a spherical 
outflow has previously been considered in the literature (e.g. 
Beloborodov 2010; Pe'er & Ryde 2011). In the limiting case 
of outflows with 9j > 5/Fo observed at 9 V <JC 0- 3 , wc obtain 
similar results. For larger viewing angles or smaller jet open- 
ing angles, geometrical broadening of the spectrum has to 
be considered. 

Although we consider a static source, as shown in §6.2 
there is a time delay associated with high-latitude photons. 
The time delay increases with decreasing observed photon 
energy. The time delay at a specific energy is longer for 
more narrowly collimated jets. For a narrow jet {9 S = 1/To) 
with typical outflow parameters characterizing GRBs, the 
time delay is « 1 s at E « 10~ 2 i5 P cak. Spectral analysis of 
prompt GRB emission using smaller time bins than this time 
delay may reveal harder spectra within the energy range 
10 _2 -E P cak < E < -Epcak than what is predicted by the static 
model. 

For this work we have considered outflows with angle 
independent luminosity, and so the shape of the Lorentz 
factor profile is fully determined by the angle dependent 
baryon loading. As shown above, the spectral slope in the 
range 10~ 2 -E poa k < E < -E pC ak is formed by photons making 
their last scattering at 8 < 5/Fo. Therefore, this assumption 
is a good approximation for jets where dL/d£l « const for 
9 < 5/Fo- This requirement is fulfilled for model JA in Zhang 
et al. (2003) close to the largest radius of the simulation 
(r = 2.1 x 10 10 cm). In our notation, r « 140, 9 S « 0.017 « 
2.5/To and p « 6.5 for model JA as estimated from Figures 
8 and 9 in Zhang et al. (2003). 

Furthermore, we consider non-dissipative fireball dy- 
namics. The dynamics of the outflow are dependent on the 
dominant form of energy carried by the jet as well as en- 
ergy dissipation. In particular, the radial scalings of the 
comoving temperature and Lorentz factor in magnetically 
dominated outflows are expected to differ significantly from 
the scalings of thermal fireballs. Giannios (2012) and Be- 
loborodov (2012) considered the scaling of -E pe ak in dissi- 
pative outflows. Heating keeps the comoving temperature 
approximately constant in the range i? p h/30 < r < 7? p h, 
in contrast to the non-dissipative outflows considered here. 
However, the framework presented in §2 is general and may 
be applied to any relativistic, optically thick outflow. 

Since the causally connected parts of the outflow are 
separated by angles « 1/F, one may consider outflows with 
Lorentz factor variations at angular scales of A9 « few/F. 
Geometrical broadening of the observed spectrum is ex- 
pected as a consequence of the beaming of photons being 
a function of angle from the jet axis, in much the same way 
as for the jets considered in this work. The spectral slope 
below the peak energy is expected to depend on both the 
typical angular scale as well as the amplitude of the Lorentz 
factor variations. 

Spectral broadening by energy dissipation in regions of 
moderate optical depth may be combined with geometri- 
cal broadening. As the observed spectrum below the peak 
energy is a superposition of comoving spectra, it is not sen- 
sitive to the exact shape of the comoving spectrum. Comp- 
tonization of the comoving spectrum by electrons which are 
heated by energy dissipation (e.g. Beloborodov 2010; Laz- 
zati & Begelman 2010) can shape the spectrum above the 
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peak energy, while geometrical broadening forms the spec- 
trum below the peak energy. 

The observed low energy photon index varies between 
bursts, forming an approximately Gaussian distribution cen- 
tered at a « — 1 with a full width at half maximum of 
« 1 (e.g. Goldstein et al. 2012). Within the framework pre- 
sented in this paper, the distribution could naturally be in- 
terpreted as a result of viewing angle variations. The ex- 
act a-distribution predicted by our model is hard to obtain, 
since it depends on detector characteristics. However, a clear 
prediction of the model is that of softening of the low en- 
ergy photon index with increasing viewing angle for jets with 
*j > few /T . 

GRB 090902B is a burst of special interest, due to its 
unusual spectral shape (Abdo et al. 2009). The prompt spec- 
trum consists of a sharply peaked Band component along 
with a wide power law. The hard low energy photon in- 
dex and the narrow width of the Band component poses 
an extreme challenge for optically thin emission models. 
Therefore, a photospheric origin of the Band component in 
this burst seems inevitable (Ryde et al. 2010, 2011; Zhang 
et al. 2011; Pe'er et al. 2012). Since the Band component 
appears close to blackbody, processes that are expected to 
modify the photospheric spectrum can be constrained for 
GRB 090902B. The rate of energy dissipation in regions of 
moderate optical depth must be relatively low as the ob- 
served Band component is not severely distortened from the 
Planck spectrum. Within the framework considered in this 
paper, further constraints can be set. The requirement for 
geometrical broadening to not significantly distort the ob- 
served spectrum in the energy range 10 _2 -E pea k < E < i? pea k 
is 9 V <JC 6 S and 6 S > 5/To. Depending on typical GRB outflow 
characteristics, the probability for such parameter combina- 
tions may be low. This would help explaining the rarity of 
similar events. 

Multiple spectral components have been clearly identi- 
fied in several GRBs after the launch of the Fermi telescope 
(Abdo et al. 2009; Ackermann et al. 2010, 2011). This was 
expected from a theoretical point of view, simply because of 
the wide energy range of the observations. Several models 
predict non-thermal spectral components as a result of ki- 
netic or magnetic energy dissipation above the photosphere 
(e.g. Meszaros et al. 2002; Zhang &: Yan 2011). If dissipation 
occurs at r ^ i? p h the thermal and non-thermal components 
may be fitted separately (such as in 090902B (Pe'er et al. 
2012)). For dissipation at r > i? ph the separation of ther- 
mal and non-thermal components may be less clear because 
of the coupling of the thermal photon field to the acceler- 
ated non-thermal electrons. In this work we consider the 
thermal component in isolation in order to demonstrate the 
geometrical broadening of the spectrum. However, a com- 
plete theory for the prompt GRB emission must explain all 
observed spectral components. 

An outstanding problem for the photospheric interpre- 
tation of prompt GRB emission is the non-thermal spectra 
commonly observed. In particular, the average low energy 
photon index is not compatible with the Rayleigh-Jeans in- 
dex (q £s — 1 as compared to a = 1). In this work we show 
that either a narrow jet (6>j < few/To) with a moderate 
Lorentz factor gradient (1 < p < 4) observed at any viewing 
angle or a wide jet observed at 6 V « 9j can naturally produce 
— 1 < Q: < —0.5 through geometrical broadening, without a 



need for additional emission processes such as synchrotron 
emission to supply photons below the thermal peak. 
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